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We study theoretically the decoherence of a gas of bosonic atoms induced by the interaction 
with a largely detuned laser beam. It is shown that for a standing laser beam decoherence coincides 
with the single-particle result. For a running laser beam many-particle effects lead to significant 
modifications. 
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In experiments with atomic Bose- Einstein condensates 
the trap which is used to isolate the atoms from 
the environment plays an important role for the physi- 
cal behaviour of the system. In a magnetic trap usually 
only atoms with a specific magnetic hyperfine quantum 
number are stored. An optical trap |^ can be used to 
confine atoms in all Zeeman sublevels. However, in order 
to preserve the atomic coherence spontaneous emission 
of photons must be avoided. A large detuning A of the 
laser beams providing the trap field reduces the excita- 
tion probability for atoms and thus reduces the number 
of spontaneously emitted photons. At the same time, a 
relatively strong intensity of the laser beams, and cor- 
respondingly a large Rabi frequency ri(x), ensures that 
the laser beams' effect is still large enough to produce a 
strong potential for the atoms. 

The aim of this paper is to examine light-induced deco- 
herence of a gas of ultracold atoms in an optical trap and 
to deduce the influence of atomic many-particle proper- 
ties on this decoherence. To do so we study solutions 
to a master equation for reduced Heisenberg operators 
describing a quantum field of bosonic two-level atoms 
interacting with a detuned laser beam and the vacuum 
fluctuations of the electromagnetic field. The technique 
of reduced Heisenberg operators has previously been em- 
ployed by Nemoto and Shibata Q and allows an elegant 
description of decoherence in a many-particle context. 
After the derivation of the master equation we will adia- 
batically eliminate the excited state to get an effective 
master equation for the atomic internal ground state. 
The main result of the paper is that many-particle ef- 
fects lead to a nonlocal modification of light-induced de- 
coherence that vanishes for laser fields with spatially ho- 
mogeneous phase (e.g., a standing laser wave) but makes 
substantial contributions for a running laser wave. 

Derivation of the master equation: We consider a 
general system that is described by the direct product of 
two Hilbert spaces TLs and TLr, where Tis describes the 
open system that we are interested in, and TLr describes 
the reservoir to which the system is coupled. It is our 
aim to find a master equation describing the evolution 
of reduced Heisenberg operators for the system only. We 
define a reduced Heisenberg operator Rs by the relation 



Rs{t) -TrRipRRit)} 



(1) 



where R{t) is the original Heisenberg operator which acts 
on Tis ® Ti-R and pR is the density matrix describing the 
state of the reservoir. It is thereby assumed that the 
(time-independent) density matrix can be written in the 
form p = PS 'S> PR so that there are initially no correla- 
tions between the system and the reservoir. The physical 
significance of the reduced Heisenberg operator is that it 
allows to describe accurately any measurement involving 
the system only, i.e., 

Trs^RipFsmP;,} = TTs{psFsRsit)Fi.} (2) 

for any two operators Fs and F'g acting on Tis only. 

To derive a master equation for Rsit) we employ 
Zwanzig's method by using that VR := Ttr{prR} 1r 
is a projection superoperator fulfihing = V. We start 
from the Heisenberg equation-of-motion 



ihdtR{t) = [R{t),H] =: CR{t) 



(3) 



where C is the Liouville superoperator. The Hamiltonian 
in Heisenberg picture is explicitely time-independent if 
it is written in terms of Schrodingcr picture operators. 
Introducing Q 1 — V and following the lines of 
Ref. one can derive the master equation for the re- 
duced Heisenberg operator. 



ihdtRsit) ^PCPRsit) 



(4) 
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dt'UQQ{t')QCVRs{t~t') 



In the derivation it was assumed that the inital 
Heisenberg operator i?(0) acts on the system only so 
that QR{Q) = 0. The superoperator Uggit') := 
exp{—it'QCQ/h} describes the unitary evolution in Q- 
subspace. 

This expression can be further simplified if the Liou- 
villean is written in the form C = Cs + Cr + Cmti where 
Cs and Cr act solely on the system and the reservoir, re- 
spectively, and /lint describes the interaction. It is then 
easy to show that QCV = QC-mtP and, if the reservoir 
is in a stationary state {[Hr,pr\ = 0), VCQ = VCmtQ- 
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Expanding Eq. (Q) to second order in >Cint we can re- 
place UQQ{t') by exp{— ii'(>C5 + Cr)/Ti\. Performing a 
Markov approximation {Rs{t~t') « Upp[t')Rs{t) under 
the integral) then results in the final form of the master 
Equation in Heisenberg picture, 

POO 

ihdtRs[t) « CsRs{t) - '-VC^tQ / dt'£int{t')Rsit) , 
^ Jo 

(5) 

with Cint{t')Rs = [Rs,Hint{t')] and Hi^tit') 
exp{-it'(£s + CR)/h}Hi^f 

Decoherence of the atomic field operator: The 
formalism of the previous section will now be applied 
to the problem of spontaneous emission in a system of 
bosonic two-level atoms. In this case the degrees-of- 
freedom of the electromagnetic field play the role of a 
reservoir to which the atoms are coupled. Since it is 
known that usually the center-of-mass motion of the 
atoms has only a small influence on the spontaneous 
emission rate we neglect it alltogether and use the 
atomic Hamiltonian 



(6) 



where loq is the resonance frequency and 'l'i(x) are the 
quantum field operators for excited {i = e) and ground- 
state {i = g) atoms, respectively. The Hamiltonian of 
the reservoir is given by 



(7) 



and the coupling between the system of two-level atoms 
and the reservoir is described by the electic-dipole cou- 
pling in rotating-wave approximation. 



E(-)(x,<')P^+'(x)e-''^°*'} , (8) 

with E*^+)(x, i') being the positive-frequency part of 
the electric field in interaction picture and P*^+)(x) := 
d5'e(x)\l/J(x) is the positive-frequency part of the polar- 
ization operator. 

Under the assumption that the initial state of the elec- 
tromagnetic field is the vacuum |0) we can insert the 
various Hamiltonians into Eq. (||) to derive 

ihdtRsit) = [Rs{t), Hs] - ihCoRsit) , (9) 
with the decoherence Liouvillean 

CdR ■■= y d3a;d3a;'|[i?,/^(-'(x)]pj+'(x')Ty(x,x') + 

p(-)(x)[i^W(x'),i?m*(x,x')}. (10) 



The quantity Tij can be calculated by standard methods 
(see, e.g., Refs. |^,^) and is given by 

/>oo 

T,,(x,x'):-/ dt'{0\El+\^,t')El-\^',Om (H) 
Jo 



diVLO 
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In the last expression we have defined r := |x — x'| and 
f := (x — x')/r. 

It is not hard to see that the imaginary part of Tij leads 
only to a (non-local) modification of the Hamiltonian Hs 
and therefore does not lead to decoherence effects. Since 
it is the latter which we are interested in we will omit the 
imaginary part. Doing so and inserting the definition of 
the polarization operator into Eq. (fol) results in 
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'{[4*i(x)*,(x)]*t(x')*e(x') + 

*t(x)f,(x)[*t(x')*e(x'),i?]}j(r,0) . (12) 



Here 70 :— |dpWg/(37r?ieoc'^) is the single-atom sponta- 
neous emission rate and the function J(r, 9) is defined 
as 



J(r,^) sm^(^)jo(^^) + 4^ (^^^^'W " l) 



(13) 



where 9 is the angle between d and f . In the context of 
master equations for reduced density matrices the deco- 
herence term (^2[) is well known |l^,|ll| and has recently 
been applied to study loading ||l2| and condensation |l^] 
of a condensate in laser light. 

Eq. ([ij) will be the starting point of our investiga- 
tion on light-induced decoherence in many-atom systems. 
Before doing so we point out that the master equations 
Eq. ([1^) and (|o|) are very similar to corresponding equa- 
tions for the reduced density matrix in Schrodinger pic- 
ture. However, there is a striking difference: For the re- 
duced density matrix ps the Liouvillean for spontaneous 
emission includes terms in which the operators P*^+^(x) 
and P(~)(x) appear to the left and to the right of the 
density matrix, respectively. This ordering is important 
as such a term describes the incoherent de-excitation of 
an atom, while the opposite ordering P(+^ (x)/5sP^^^ (x) 
would describe an incoherent excitation of an atom. But 
it is exactly the latter type of ordering which appears in 
Eq. (px|). The reason is that Eq. ( p^ ) does not describe 
the evolution of the density matrix in Schrodinger picture 
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but that of an operator in the Heisenberg picture. When 
appUed to a Heisenberg operator the "wrong" ordering 
gives indeed the correct time evolution. 

A gas of atoms in largely detuned laser light: 
To study the decoherence of atoms in a largely detuned 
laser beam we adiabatically eliminate the excited state 
by means of (see, e.g., Ref. Q) 



(14) 



Inserting this into Eq. (12) leads to a Liouvillean 
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3 3 , r!*(x)»(xO 
a xd X — 



J(r, 6) X 



{[i?,p<,(x)]/5<,(x') -/5<,(x)[i?,p,(x')]} (15) 

with pg(x) :=*t(x)*g(x). 

Eq. ([1^) describes the decoherence of any Heisenberg 
operator acting on the ground state of the atoms in the 
presence of a laser beam. Its physical interpretation is 
that the laser beam excites an atom at position x' to 
a virtual state. The atom then spontaneously emits a 
photon which is absorbed by another atom at position 
X. Finally, this atom makes an induced emission of a 
photon into the mode of the laser beam. A simple ap- 



plication of Eq. (15) is the demonstration of lack of de- 
coherence for the number density operator /5g(x). Since 
[pg(x), /5g(x')] = one can immediately infer that 



(16) 



i.e., there will be no decoherence of the number density 
at a point x. This is a reasonable result since sponta- 
neous emission should not change the population density 
if the atomic center-of-mass motion is neglected. 

To study the phase decoherence we consider the time 
evolution of the field operator ^g(x) itself. Introducing 
the operator 



Q(x) 
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dVj(r,6l)17(x')f^*(x)pg(x') (17) 



one can show that (£g)"\t'g(x) = (£Q)"*I'g(x) with 

CqR := {i^)R - RQ{x) . (18) 

We therefore can replace Cg by Cq for the special case 
of the field operator. The solution of the equation of mo- 
tion dfR = —CqR with the initial condition R{t = 0) = 
^g(x) then can be written as 



i?(t) =e*'3^W|'g(x)e-*'3W 



(19) 



This expression can be further reduced by writing it in 
the form R{t) = exp[-t((5-gt)] exp[tQ]*g(x) exp[-t(5]. 
The first exponential corresponds to a unitary opera- 
tor. The other two exponentials represent a non-untitary 



transformation of the field operator. It is easy to prove 
jlSf that this non-unitary transformation can be reduced 
to 



=*Q(x) 



*s(x)e 



-tQ(x) _ 



exp 



lat |r!(x)|^ 



*s(x) (20) 



This result just corresponds to the single-particle 
decoherence-effect of spontaneous emission: at point x 
the probability for the atoms to become excited is, within 
the adiabatic approximation, given by |il(x)p/A2. The 
excited atoms then do spontaneously (i.e., incoherently) 
decay at a rate 70. The rate 7o|il(x)p/A2 therefore de- 
termines the time scale of light-induced decoherence at 
point X. 

It is surprising that the transformation ( pO| ) docs not 
describe the full decoherence of the field operator but 
that one also has to take into account the unitary op- 
erator Uq := exp[— i((3 — Q^)]. If Q is Hermitean the 
unitary operator becomes trivial and the total effect of 
decoherence reduces to that described by Eq. (po|). This 
happens whenever the phase of the laser beam is homo- 
geneous, for instance in a standing-wave laser field with 
fixed polarization. The reason is the cancellation between 
the contributions of counterpropagating light waves. In 
a running wave, however, Q is not Hermitean and many- 
particle effects do appear. 

To analyse the effect of the unitary operator Uq it 
is convenient to consider a fully coherent ensemble of 
atoms, e.g., a Bose-Einstein condensate at temperature 
T = 0. In the mean-field description the condensate can 
be described (see, e.g., Ref. [||) by a coherent state |q;) = 
exp(-|ap/2) exp(Q!aJ)|0), where = / d^x(^o(x)5'j,(x) 
is the creation operator for atoms in the condensate mode 
</?o(x). The number |ap is the number of atoms in the 
condensate. Our task is then to calculate the expectation 
value {a\R{t)\a) which, after some elementary algebra, 
can be cast into the form 



{a\R{t)\a) ~ Q!(^o(x)exp 



7oi |r!(x)|2 



|0) 



(21) 



where a\{t) is defined as Uqcl^Uq. With techniques 
similar to that used to derive Eq. ( [20| ) one can prove 
that a\{t) is in fact the creation operator for another 
mode, aj(t) = / d"^a;(^o(x, t)5'J(x), whose mode function 
is given by 



(^o(x',i) = (/?o(x')exp 



- -170^ J(x - x') X (22) 
r2(x')r2*(x) - f7*(x')f^(x) 
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Inserting this expression into Eq. (|2^) leads to the final 
result for the light-induced decoherence of a condensate 
at position x, 
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{a\R{t)\a) — Q;(po(x)exp 



lot |r!(x)|^ 



A2 



(23) 



exp 



d3xV*(x')^o(x',<)-l^ 



Eq. ( P3|) is the main result of this paper and describes 
the total effect of light-induced decoherence on a Bose- 
Einstein condensate of two-level atoms. It demonstrates 
that there are non-local contributions to the atomic de- 
coherence. To estimate the order-of-magnitude of these 
non-local terms we assume for simplicity a Gaussian 
condensate wave function of width w, i.e., (po{^) = 
expl—x^ / {2w'^)]/ (w'^^'^TT^^^) . In addition we consider a 
running-wave laser field of the form ^(x) = &ilo exp(ifcz). 
The analysis of Eq. ( p3| ) is then reduced to a numerical 
evaluation of a complicated integral. To gain more physi- 
cal insight we expand the exponent appearing in Eq. ( p^ ) 
to second order in t. Introducing the decoherence rate 
7£) := 7o|rio|'^/(2A2) we then can write the result ( |2^ ) 
as 

{a\R{t)\a) « ck(^o(x) exp — 7d^(1 + 3iiVAA(z, w, k^)) 



(24) 



Here N\ :— \a(po{0)\'^XQ is the number of condensed 
atoms in a cube of dimension Aq at the center of the 
condensate. Aq is the wavelength of resonant light. The 
functions A{z, w, kL) and B(z^ w, k^) depend on the po- 
sition X = zGz at which the decoherence is observed, the 
width w of the condensate, and on the wavevector /cl of 
the running laser wave. For w — lOOc/cjo ^md a detuning 
of A = 1 GHz they are shown in Figs. and ||, respec- 
tively. A only has a strong dependence on k^ for extreme 
detunings above 1 THz while B is rather insensitive to 
the detuning. 

Eq. ( p^ ) demonstrates that the additional unitary op- 
erator Uq does not lead to an additional exponential de- 
cay of the coherences since to first order its effect is only 
a position-dependent phase shift which is proportional to 
A. For longer times {t w 1/jd) the many-particle deco- 
herence effect can become even larger than the single- 
particle contribution. Since in the example given above 
the function B is roughly proportional to the conden- 
sate density one arrives at the convincing result that the 
many-particle decoherence grows with the atomic den- 
sity 

In conclusion, we have shown that the decoherence of 
a gas of ultracold atoms in a detuned laser beams is spa- 
tially varying and depends on whether the phase of the 
laser beam is spatially homogeneous. If this is not the 
case the one-particle spontaneous emission process will 
be accompagnied by a significant many-particle contribu- 
tion with a non-exponential decay rate. This additional 
effect corresponds to the expectation value of the unitary 
operator Uq. 
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One has to show that the left-hand and the right-hand 
side of Eq. (|^) are the same at t = and that they fulfill 
the same diflterential equation. 




FIG. 1. The function A{z,w,kL), which is related to a 
light-induced phase shift in the condensate, as a function of 
the position z along the z-axis. The width w is 100/(27r) times 
the resonant wavelength, and the detuning was chosen to be 
A = IGHz. 
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FIG. 2. The function B{z,w,kL), which is related to 
the many-particle decoherence effect, for the same pa- 
rameters as in Fig. 0. 



